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Abstract. In recent years, we used lattice QCD to calculate some quantities that were 
unknown or poorly known. They are the q 2 dependence of the form factor in semileptonic 
D — > Klu decay, the leptonic decay constants of the D + and D a mesons, and the mass of the 
Be meson. In this paper, we summarize these calculations, with emphasis on their (subsequent) 
confirmation by measurements in e + e~ , -yp and pp collisions 



1. Introduction and Background 

The central theme of elementary particle physics is to find new interactions of matter, energy, 
space and time. When the matter in question is the quarks, one is faced with quark confinement: 
quarks never appear freely; they are always bound inside hadrons — baryons like the proton, or 
mesons like the pion or kaon. In the Standard Model of elementary particles, confinement is a 
phenomenon of quantum chromodynamics (QCD), the gauge theory of the strong force. 

Since only hadrons can be detected, the effects of quark confinement must be calculated with 
QCD, before the experimental data can be interpreted in terms of quarks. In many cases, the 
best technique for doing the calculations is to formulate QCD on a space-time lattice. Lattice 
QCD and the Feynman path integral reduce the problem of to an integral whose dimension 
scales as ./V 4 , where N is the (linear) lattice size. The problem cries out for supercomputing. 

In recent years, lattice QCD has reached the stage where many calculations of hadron masses, 
mass splittings and operator matrix elements agree with experimental measurements. The key 
has been the inclusion of sea quarks, which are pairs of virtual quarks swirling around inside 
hadrons. The progress has been especially striking when the sea quarks are implemented as 
staggered quarks, using an action designed to reduce discretization effects. 

One ingredient of these calculations is controversial. Sea quarks are always computationally 
demanding, although staggered quarks are by far the fastest. Staggered quarks introduce some 
extra unwanted quarks, however. The computer algorithms [2] and subsequent analysis of the 
numerical data j3J remove them, but do so differently for valence and sea quarks. The difference 
can lead to violations of unitarity. In the cases discussed here, it is plausible that such effects 
are small, but a proof is not yet at hand [I]. 

Less controversial is the treatment of heavy quarks. In practice, the lattice spacing is not 
small enough to resolve the Compton wavelength of charmed and b quarks. Fortunately, 
chromodynamics at this length scale is simple enough to factor it out from the computer 
simulation, and several methods exist [Sj. Nevertheless, it is good to have check. 



In this paper, we discuss three topics: the normalization and g 2 -dependence of the D — » Klv 
form factor; the decay constants of the D + and D s mesons; and the mass of the B c meson. Each 
of these lattice-QCD calculations was subsequently confirmed by experimental measurements, 
satisfying a long-standing demand of experimental physicists [§]. The quantities discussed here 
were ideal candidates: they are straightforward to compute; they test the controversial aspects 
in complementary ways; and the first "good" experimental measurements were expected on the 
same time scale. The success of the predictions is extremely encouraging. In particular, the 
calculations for D mesons are, in lattice QCD, similar to those for B mesons, whose b quarks 
are considered likely to exhibit new, non-Standard interactions. 

2. Semileptonic D Decays 

Semileptonic decays such as D — > Klv proceed as follows. A quark (in this case, a charmed 
quark) emits a virtual W boson, thereby turning into a quark of a different flavor (in this case, 
a strange quark). The W immediately disintegrates into a lepton-neutrino (lu) pair. The rate 
depends on q 2 , which is the invariant-mass-squared of Iv. Some of the q 2 dependence stems from 
QCD through a function called a form factor (in this case, denoted /+(g 2 )). The momentum 
transfer q 2 falls in the range < q 2 < q 2 nax = (mo — tuk) 2 - In lattice QCD, discretization effects 
are smallest when the spatial momentum p of the kaon is small, which puts q 2 close to g^ ax . 

Experiments usually measure the branching fraction and quote the normalization /+(0), 
after making assumptions about the q 2 dependence. While our results were still preliminary 0, 
experimental results came out for the normalization of D — > Klv [S] and D — » ttIv The 
agreement with our final results ^01 is excellent. For example, we find /+ ~ ^(0) = 0.73(3)(7) JOJ 
while the BES Collaboration measures f+^ K (0) = 0.78(5) 0. 

In principle, the shape of the form factors can be computed directly in lattice QCD. In 
practice, we calculated at a few values of p and used a fit to the Ansatz of Becirevic-Kaidalov 
(BK) 11 to fix the q 2 dependence. It was important, therefore, to measure the q 2 dependence 
experimentally. In photoproduction of charm off fixed nuclear targets, the FOCUS Collaboration 
was able to collect high enough statistics to trace out the q 2 distribution of the decay [T2] . 
This setup does not yield an absolutely normalized branching ratio, so one is left to compare 
U(q 2 )/U(0). 

In Fig. ^a) we plot our result for f + (q 2 )/f + (0) vs. q 2 /m 2 D „. The errors from /+(0) must 
be propagated to non-zero q 2 , so for f + (q 2 )/f + (0) the errors grow with q 2 . Figure ^ shows 1-a 




Figure 1. Form factor for D — > Klv vs. q 2 /m 2 D *\ (a) shape f+(q 2 )/ /+(0) compared with 
FOCUS ^2]; (b) shape and normalization f+{q 2 ) compared with Belle 14 . 



bands of statistical (orange) and all uncertainties (yellow) added in quadrature |13) . As one 
can see, the q 2 dependence of lattice QCD (curve and error band) and data from the FOCUS 
experiment (points) agree excellently, although the uncertainties are still several per cent. The 
FOCUS results appeared two months after the lattice calculation. More recently, the Belle 
Collaboration at the e + e~ collider KEK-i? measured the shape and normalization of the form 
factor in a single experiment ^3]. In Fig. E^b) we compare our result for f+(q 2 ) with Belle. The 
color code for the lattice QCD error bands is as before, and now depict q 2 dependence of the 
lattice-QCD errors in a realistic way. 



3. Leptonic D Decays 

We also considered the leptonic decay of charmed mesons, D + — > Iv and D s — > Iv. Here the 
quark and antiquark in the meson merge into a virtual W, which disintegrates into Iv. The 
QCD influence is a single number (for each meson), called decay constants and denoted fjj+ 
or fj} s . At Lattice 2004 |15j . we presented preliminary results for f D + and fjj s , based on one 
lattice spacing, a « 0.125 fm. Our extended the running to two other lattice spacings. Details 
are given in the ensuing publication JTCj. We find 

f D + = 201±3±17MeV, (1) 
f D „ = 249 ± 3 ± 16 MeV, (2) 

where the first error is from finite Monte Carlo statistics, the second is a sum in quadrature 
of several systematics. A conservative (but not naive) estimate of heavy-quark discretizations 
effects is the second largest (largest) systematic on f D + (/d s ). 

Figure [3 shows the rif dependence of the decay constants. Quenched (rif = 0) results vary 
widely, but we show one J3| carried out with similar choices for other aspects of the calculations. 
One sees a trend of fn 3 to increase with rif. A similar comparison of /d+, in Fig. 12(b), is less 
instructive, but shown for completeness. 

The CLEO-c Collaboration ^2] and the BaBar Collaboration [20] have measured 

f D + = 223 ±17 ±3 MeV CLEO-c PU, (3) 
f Ds = 279 ±17 ±20 MeV BaBar (20], (4) 

respectivelty at the CESR and PEP-II e + e~ colliders. At the 1-a level, the agreement with 
lattice QCD is fine. Even more compelling is the ratio R^u = /d s : 

R d/s = 0.786 ± 0.042 lattice QCD, (5) 
= 0.779 ±0.093 CLEO/BaBar, (6) 
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Figure 2. Dependence of (a) fjj a and (b) f D + on the number rif of sea flavors. Quenched (rif = 
0) rif = 2 ^H]; rif = 3 JS]- Solid (dashed) error bars are statistical (statistical±systematic). 



in which several uncertainties from lattice QCD cancel. Experimental and lattice-QCD 
uncertainties are reducible, so the test will sharpen over the coming few years. 



4. Mass of the B c Meson 

The pseudoscalar B+ meson is the lowest-lying bound state of a charmed quark and a b quark. 
The CDF Collaboration [2J first observed it during Run I of the Tevatron pp collider in the 
semileptonic decay B£ — > J/tpl + i^. But it was clear that Tevatron Run II detectors would 
be able to reconstruct hadronic modes, such as B+ — > J/tfjTT + , which give much much better 
precision on m,B c 22 j. At Lattice 2004 we presented results in nearly final form [20], and posted 
the final results on the arXiv in mid-November |24j : 

m Bc = 6304 ± 12± 1 g MeV, (7) 

where the last error is a rough estimate of residual heavy-quark discretization effects. Soon 
afterwards, CDF announced their precise mass measurement [25| : 

m Bc = 6287 ± 5 MeV, (8) 

which agrees with Eq. (J7J) at slightly more than 1-a. 

Two comments are in order. First, the agreement at the gross level of the calculation with 
experiment shows that discretization effects are well under control with the heavy-quark methods 
of choice. These are lattice NRQCD [27] and the Fermilab method [25 , which are based on 
effective field theories for heavy quarks [291 I30j . Indeed, as seen in Fig. |3^a), almost no lattice 
spacing dependence is seen in the splitting A^y = "1b c — {m^ + my)/ 2 that is at the crux 
of the calculation [20]- Moreover, it is striking how little the splitting A^,x changes when sea 
quarks are included. FigureCUb) compares Eq. |j7J) with an old quenched calculation [20] ( an d the 
measurement [221 )• The solid error bar shows the non-quenching errors, and the dashed includes 
the estimate of the quenching error. The inclusion of sea quarks has reduced the splitting by a 
factor of three or four, bringing an essentially discrepant result into agreement. 



5. Conclusions 

In the past year, several lattice-QCD calculations have been confirmed by experiment. FOCUS 
|12j and Belle Jl] confirmed the g 2 -dependence of the D — > Klv form factor [10] : CLEO- 
c ^3] and BaBar [20] respectively confirmed the D + and D s decay constants |16| : and 
CDF 25 confirmed the mass of the B c meson [24] . To obtain these results it is essential 
to have heavy-quark discretization effects under control, as one expects from theoretical 
foundations (27J [2H1 1201 EO] • Furthermore, the comparison of quenched QCD, QCD with 2+1 
staggered flavors, and experiment shows that sea quarks are needed to obtain agreement, and 
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Figure 3. (a) Dependence of the splitting A^y on the lattice spacing a. (b) Comparison of the 
quenched f2B., nf = 2 + 1 [21] and experimental [20] values of ms c ; the dashed line denotes the 
baseline (m^p + my)/ 2. 



that staggered quarks (in these cases) capture the needed effect. The results are promising for 
the search for new b quark interactions, because a straighforward change to the D form factors 
and decay constants yield the corresponding results for B mesons. These are a key element to 
enable the experimental search for new phenomena in quark- flavor physics 6 . 
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